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ABSTRACT 

A recent re-determination of the non-thermal component of the hard X-ray to soft 7-ray emission from 
the Galactic ridge, using the SPI instrument on the INTErnational Gamma-Ray Astrophysics Labo- 
ratory (INTEGRAL) Observatory, is shown to be well reproduced as inverse-Compton emission from 
the interstellar medium. Both cosmic-ray primary electrons and secondary electrons and positrons 
contribute to the emission. The prediction uses the GALPROP model and includes a new calculation 
of the interstellar radiation field. This may solve a long-standing mystery of the origin of this emission, 
and potentially opens a new window on Galactic cosmic rays. 
Subject headings: elementary particles — cosmic rays — gamma-rays: theory 



1. INTRODUCTION 

t^J- ' The Galactic ridge is known to be an intense source 
[--"-. of continuum hard X- and 7-ray em ission. The hard X - 
t — ■ ray emission was discovered in 1972 (Blea ch et al.lll972j) . 
t-H ' and interstellar emission has subsequently been observed 
^-J by HE AO-1, Tenma (ASTRO-B) ASCA, Ginga, RXTE 
q ■ OSSE (iKinzer et all Il999l [2 001; K oyama et all [l98l 
OO ! iPurcell et al.l 11996: W orrall et al.l 1 19821 ). and most re- 
f^i , cently by Chandra and XMM-Newton. The 7-ray ob- 
servations started with the OSO-III satellite in 1968, fol- 
. £ ' lowed by SAS-2 in 1972, COS-B (1975-1982) and COMP- 
TEL and EGRET on the CGRO (1991-2000). With 
COMPTEL and EGRET the improvement in data qual- 
ity was sufficient to allow such studies to be performed 
in much greater detail. The Galactic diffuse emission is 
a major study objective for INTEGRAL (in orbit since 
2002) and the GLAST LAT (to be launched in 2008) 
(|Michelsonll2007t lRitzll2007f ). Each of these experiments 
represents a significant leap forward with respect to its 
predecessor. 

Continuum emission of diffuse, interstellar nature is ex- 
pected in the hard X-ray and 7-ray regime from the phys- 
ical processes of positron annihilation (through interme- 
diate formation of positronium), inverse-Compton (IC) 
scattering and bremsstrahlung from cosmic-ray (CR) 
electrons and positrons, and via decay of neutral pions 
produced by interactions of CR nuclei with the inter- 
stellar gas. Positron annihilation in flight (continuum) 
may contribute in the few MeV range (|Beacom fe; Yiiksel 
2006). For the non-positronium continuum, hard X- 

1 Also Kavli Institute for Particle Astrophysics and Cosmology, 
Stanford University, Stanford, CA 94309 



rays from bremsstrahlung emission imply a luminosity 
in CR electrons w hich is unacceptably large (see e.g. 
iDogiel et al.ll2002"al ). Composite models have been pro- 
posed which incorporate thermal and nonthermal com- 
ponents from electrons accelerated i n supernovae or the 
ambient interstellar turbulence (Valinia et al. 2000b). At 
MeV energies the o rigin of the emission is also uncertain 
([Strong et all l2TXJ0h . 

Alternatively, the origin of the ridge emission could 
be attributed to a population of sources too weak to 
be detected individually, and hence would not be truly 
interstellar. In general, 7-ray telescopes have inade- 
quate spatial resolution to clarify this issue. For X- 
rays, important prog ress in this a rea w as made by ASCA 
(|Kaneda et al.lll997f ) and Ginga (jYamasaki et al.l ll997). 
More recently, high-res olution imaging in X-rays (2- 
10 keV) with Chandra ([Ebisawa et al.ll200lL 120061 ) has 
claimed to prove the existence of a truly diffuse compo- 
nent. Similiarly, it h as been claim ed from an analysis of 
XMM-Newton data ([Hands et al.ll2004ft that 80% of the 
Galactic-ridge X-ray emission is probably diffuse, and 
only 9% can be accounted for by Galactic sources, the 
rest b eing extragalactic in na ture. However, more re- 
cent ly IRevnivtsev et all (I2006D . with RXTE PCA data, 
and lKrivonos et all (|2007] )7using INTEGRAL/IBIS and 
RXTE data, argue convincingly that below 50 keV all 
the "diffuse" emission can be accounted for by a Galac- 
tic population of so urces, mainly magnetic catacl ysmic 
variables; see al s o IRevnivtsev fc Sazonovl (|2007l ) and 
IRevnivtsev et all (|2007f) . 

At higher energies, an extensive study of the diffuse 
Galactic 7-ray emission in the cont ext of CR propaga- 
tion models has been carried out bv lStrong et alj (|2000l 
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l2004bh . This study confirmed that models based on 
locally measured electron and nuclei spectra and syn- 
chrotron constraints are consistent with 7-ray measure- 
ments in the 30 MeV - 500 MeV range; outside this 
range deviations from the data are apparent. The puz- 
zling excess in the EGRET diffuse emission data above 
1 GeV relat i ve to that expected (jHunter et al.l 119971 : 
IStrong et al.l [2000) has shown up in all models that 
arc tuned to be consistent wit h the locally measured 
CR nuclei and elect ron spectra (Moskalc nko et all 12004 
IStrong et al.ll2004bh . The excess has shown up in all di- 
rections, not only in the Galactic plane. This implies that 
the GeV excess is not a feature restricted to the Galactic 
ridge or the gas-related emission. A simple re-scaling of 
the components (7r°-decay, IC, bremsstrahlung) does not 
improve the fit in any region, since the observed peak 
is at an e nergy higher than the 7r ° -peak. For re c ent re - 
views see lMoskalenko et al.l (|2004D . IStrong et al.l (|2007h . 
and references therein. 

Assuming that the GeV excess is not an instrumen- 
tal artefact 2 , the so-called "optimised model," which ex- 
plains the GeV excess in terms of CR in tensity varia- 
tions in the Galaxy, has been proposed by IStrong et al.l 
(2004b). It reproduces the spectrum of the diffuse 7-rays 
in all directions, as well as the latitude and longitude 
profiles for the whole EGRET energy range 30 MeV - 50 
GeV at the cost of relaxation of the restrictions imposed 
by the measurements of local CR proton and electron 
spectra. At lower energies, the predictions of this model 
have never been tested because of the lack of good data. 

The study of the Galactic-ridge continuum X-ray emis- 
sion is a key goal of the INTEGRAL mission. The 
high spectral resolution combined with its imaging ca- 
pabilities promises new insights into the nature of this 
enigmatic radiation. Previous work based on initial, 
smaller sets of INTEGRAL/SPI observations have re- 
ported the detection of diffuse emis sion at a level con- 
sistent with previous experiments (jStrong et al.l 120031 : 
IStrong|[2003h . However, statistical and systematic errors 
were large, due in part to the uncertainty in the point- 
source contribution. Meanwhile, a new analysis of IN- 
TEG RAL/IBIS data (|Lebrun et all 12004 iTerrier et all 
l2004h showed that, up to 100 keV, indeed a large frac- 
tion of the t otal emission from th e inner Galaxy is due 
to sources. IStrong et al.l (|2004bl ) used the source cat- 
alogue from this work (containing 91 sources) as input 
to SPI model fitting, giving a more solid basis for the 
contribution of point sources in such an analysis. This 
exploited the complementarity of the instruments on IN- 
TEGRAL for the first time i n the context of diffuse emis- 
sion. The SPI analysis by iBouchet et alJ (|2005l ) gave 
a rather lower 50- 1000 keV power-law continuum than 
IStrong et al.l ((2005), but the er rors were large in the early 
datasets. Now, a new analysis (jBouchet et al.ll2008h with 
3 times as much SPI exposure gives better statistics, 
background handling, and point-source subtraction. 

In the present paper we focus on energies above 50 keV 
where sources do not appear to be important because of 
the rapid cut-off in the spectra of the majority, and the 
relatively small number of hard-spectrum sources. We 

2 A discussion of uncertainties and possible sources of error asso- 
ciated with determining the diffuse Galactic 7-ray e mission using 
EGRET data is available in Moskalcnko et al 



use the GALPROP model together with a new model for 
the Galactic interstellar radiation field (ISRF) to solve a 
long-standing mystery of the origin of the hard X-ray 
emission - IC emission from CR electrons and positrons 
- and to build a model of the Galactic diffuse emission in 
the energy range from keV to TeV energies, thus covering 
more that 10 orders of magnitude in energy. 

Primary CR electrons are directly accelerated in CR 
sources like supernova remnants or pulsars. Secondary 
electrons and positrons are produced via interactions of 
energetic nuclei with interstellar gas, and are usually con- 
sidered a minor component of CRs. This is indeed the 
case in the heliosphere where the positron to all-electron 
ratio is small at all energies, e + /(e+ + e~) to t ~ 0.1. 
However, as we will show, the combined secondary elec- 
tron/positron flux in the interstellar medium (ISM) is 
more than half of the primary CR electron flux at ~1 
GeV energies and below. This leads to a considerable 
contribution by secondary positrons and electrons to the 
diffuse 7-ray flux via IC scattering and bremsstrahlung 
and significantly, by up to a factor of ~2, increases the 
flux of diffuse Galactic emission below ^100 MeV. Sec- 
ondary positrons and electrons are, therefore, directly 
seen in hard X-rays and 7-rays. 

2. GALPROP CODE 



The GALPROP code ((Strong fc Moskalenkdl998h was 
created to enable simultaneous predictions of all rel- 
evant observations including CR nuclei, electrons and 
positrons, 7-rays and synchrotron radiation. 

We give a very brief summary of GALPROP; for 
details we refer the reader to th e relevant papers 
(Moskalcnko & Strong 1998, 2000; Moskalcnko et al. 
20021: iPtuskin et all 120061 : IStrong fc Mos kalcnko 1998; 
Strong et aiT TToOO. 2004t]) and the dedicated website 3 . 
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The GALPROP code solves the CR transport equation 
with a given source distribution and boundary condi- 
tions for all CR species. This includes a galactic wind 
(convection), diffusive reacceleration in the ISM, energy 
losses, nuclear fragmentation, radioactive decay, and pro- 
duction of secondary particles and isotopes. The nu- 
merical solution of the transport equ ation is based on 
a Crank-Nicholson ((Press et al.l I1992T ) implicit second- 
order scheme. The spatial boundary conditions assume 
free particle escape. Since the grid involves a 3D (R, z,p) 
or 4D (x,y,z,p) problem (spatial variables plus momen- 
tum) "operator splitting" is used to handle the implicit 
solution. For a given halo size the diffusion coefficient, 
as a function of momentum and the reacceleration or 
convection parameters, is determined by the boron-to- 
carbon ratio data. If reacceleration is included, the 
momentum-space diffusion coefficient D pp is related to 
the spatial coefficient D X x i— 0DqP S ) ((Berezinskii et all 
[19901 : iSeo fc Ptuskinl Hool . where 6 = 1/3 for a Kol- 
mogorov spectrum of interstellar turbulence or S = 1/2 
for a Kraichnan cascade, p is the magnetic rigidity, 
Dq is a constant, and (3 = v/c. Production of sec- 
ondary positrons and electrons is calculated using a for- 
malism described in Moskalcnk o fc Strong! ((1998 ) with 
a correction by iKelner et al.l (J2006). The 7-rays are 
calculated using the propagated CR distributions, in- 
cluding a contribution from secondary particles such 

' 5 http://galprop.stanford.edu 
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Fig. 1. — Spectral energy distribution of the MW ISRF in the Galactic plane. Line colouring: black, R = kpc; blue, R = 4 kpc; red, 
R = 8 kpc; magenta, R = 12 kpc. Left: maximum metallicity gradient; Right:, minimum metallicity gradient. The cosmic microwave 
background (CMB) is included in both figures and dominates the SED for wavelengths A > 600 /an. 



as positrons and electrons from inelastic processes in 
the I SM that inc r eases the 7-ray flux at MeV ener- 
gies (jStrong et al.l l2004br ). Gas-related 7-ray intensi- 
ties are computed from the emissivities as a function of 
(R, z, E~f) using the column densities of H 1 and H2 for 
galactocentric annuli based on 21-cm and CO surveys 
included in the GALPROP model. Neutral p ion pro- 
duction is calculated usin g the m ethod given by iDermerl 
(|1986allbh as described in IMoskalenko k Stron g (1998), 
or usi ng a parameterisation developed by iKamae et al.l 
(l2005f): bremsstra h lung is calculated using a formalism 
bv iKoch fc Motzl (|1959f ) as described in I Strong et all 
(2000). The IC scattering is treated using the appro- 
priate cross section for a n anisotropic radiation field de- 
veloped by IMoskalenk o k Strong! (j2000f ) using the full 
angular distribution of the ISRF. 

Cross-sections are based on the extensive LANL 
database, nuclear codes, and parameterisations 
(|Mashnik et al.l l2004f ) . Starting with the heaviest 
primary nucleus considered (e.g. Ni) the propagation 
solution is used to compute the source term for its 
spallation products, which are then propagated in turn, 
and so on down to protons, secondary electrons and 
positrons, and antiprotons. The inelastically scat- 
tered protons and antiprotons are treated as separate 
components (secondary protons, tertiary antiprotons). 
In this way secondaries, tertiaries, etc., are included. 
(Production of 10 B via the 10 Be-decay channel is 
important and requires a second iteration of this pro- 
cedure.) GALPROP includes K-capture and electron 
stripping processes, where a nucleus with an electron 
(H-like) is considered a separate species because of the 
difference in lifetime, and knock-on electrons. Primary 
electrons are treated separately. Normalisation of 
protons, alphas, and electrons to experimental data 
is provided (all other isotopes are determined by the 
source composition and propagation). Gamma rays 
are computed using interstellar gas data (for 7r°-decay 
and bremsstrahlung) and the ISRF model (for IC). The 
synchrotron emission is computed using the Galactic 
magnetic field model. Spectra of all species on the 
chosen grid and the 7-ray and synchrotron skymaps are 



output in standard astronomical formats for comparison 
with data. Recent extens ions to GALPROP include 
non-linear wave damping (jPtuskin et alJ l2006h and a 
dark matter package to allow for the propagation of 
WIMP annihilation products and calculation of the 
corresponding synchrotron and 7-ray skymaps; an 
interface between GA LPROP and the DarkSUSY code 
(|Gondolo et alJ I2004J ) will be implemented in the near 
future to allow direct calls of GALPROP from within 
DarkSUSY. 

The optimised model ([Strong et al.ll2004bh is used to 
calculate the diffuse emission in the range 10 keV - TeV 
energies. The CR source distribution is b ased on the 
Galactic pulsar distribution (Lorimer 2004), while the 
Yco- iac tors, Yco = N(R 2 )/Wco, are variable, increas- 
ing towards the outer Galaxy, and fully compatible with 
the expected variati ons based on the m etallicity gradi- 
ent and COBE data (| Strong et al.ll2004cf ) . Such a model 
reproduces the diffuse Galactic 7-ray emission for the 
whole sky as well as the radial gradient of diffuse Galac- 
tic 7-ray emissivity. 

3. INTERSTELLAR RADIATION FIELD 

The Galactic ISRF is the result of emission by stars, 
and the scattering, absorption, and re-emission of ab- 
sorbed starlight by dust i n the ISM. The m ost de- 
tailed calculation to date (St rong et al.l l2000f ). which 
includes spatial and wavelength dependenc e over the 
whole Galaxy, has been widely used. The I Strong et al.l 
(2000) model uses emissivities based on stellar popu- 
lation s based on COBE/DI RBE fits bv iFreudenreichl 
(|1998h and the SKY model of lWainscoat et al.l (|1992l ) to- 



gether with C OBE /DIRB E derived infrared emissivities 
(|Dwek et all 119971; [Sodroski et al.lll997l ). Subsequent to 
this work new relevant astronomical information on stel- 
lar populations, Galactic structure, and interstellar dust 
has become available, motivating a re-evaluation of the 
ISRF. We briefly describe our calculatio n of the ISRF ; 
furt her details can be f ound in IMoskalen ko et alJ (2006) 
and lPorter et alJ (|2006f ). 

The fundamental factors influencing the ISRF are the 
luminosity distribution from the stellar populations of 
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Fig. 2. — Spectra of CR electrons and positrons in the Galactic plane, as predicted by the adopted optimised GALPROP model. Left: 
Total (primary + secondary) and secondary electrons; Right: Secondary positrons. Interstellar spectra (IS): R = kpc (red long dashes), 
R = 4 kpc (blue short dashes), R = 8.5 kpc (black solid), also shown modulated to 600 MV. Secondary electrons are shown separately as 
magenta lines (IS and modulated) on the left panel at R = 8.5 kpc. 



the Galaxy and the radiative transport of the star light 
through the ISM. The interstellar dust absorbs and scat- 
ters the star light in the ultraviolet (UV) and optical, 
and re-emits the absorbed radiation in the infrared. 

In our model, we represent the stellar distribu- 
tion by four spatial components: the thin and thick 
disc, the bulge, and t he sphe roidal halo. We follow 
iGarwood fe Jonesl (|1987f ) and iWainscoat et all (|1992f ) 
and use a table of stellar spectral types comprising nor- 
mal stars and exotics to represent the luminosity function 
(LF) for each of the spatial components. The spectral 
templates for each s tellar type are taken from the semi- 
empirical library of iPic klcs (1998J). The normalisations 
per stellar type are obtained by adjusting the space densi- 
ties to reproduce the observed LFs in the V- and K-band 
for the thin disc. The LFs for the other spatial compo- 
nents are obtained by adjusting weights per component 
for each of the stellar types relative to the normalisations 
obtained for the thin disc LF. 

We assume a dust model including graphite, poly- 
cyclic aromatic hydrocarbons (PAHs), and silicate. Dust 
grains in the model are spherical and the absorption 
and scattering efficiencies for graphite , PAH s, and sili- 
cate grains are taken from lLi fe Drainel (|2001[ ). The dust 
grain abundance and s ize di stribution are taken from 
IWeingartner fc Drainel (|2001h (their best fit Galactic 
model). We assume a purely neutral ISM. We consider 
only coherent scattering, and a Henyey-Greenstein an gu- 
lar distribution function (|Henvev fc Greensteinlll94ll ) is 
used in the scattering calculation. The stochastic heat- 
ing of grains smaller than ~0.1 p,m is treated u si ng th e 
"thermal continuous" approach of iDraine fc Lil (|2001f h 
we calculate the equilibrium heating of larger dust grains 
by balancing absorption with re-emission as described by 
iLifc Drainel (|2001l) . 

Dust follows the Galactic gas distribution and we 
assume uniform m ixing between the two in the ISM 
(|Bohlin et al.l ll978). The dust-to-gas ratio scales with 



the Galactic metallicity gradient. Estimates for the 
Galactic [O /H] gradient vary i n the range 0.04 — 0.07 
dex kpc - (|Strong et al.ll2004cl and references therein). 
The variation of the metallicity gradient influences the 
redistribution of the mainly UV and blue component of 
the ISRF into the infrared: increased metallicity implies 
more dust, which enhances the absorption of the star 
light. The variation in the infrared component affects 
the emission in the hard X-rays (see below). Therefore, 
we consider two ISRFs corresponding to a maximal case 
of 0.07 dex kpc -1 and a minimal case with no gradient. 

The ISRF is calculated for a cylindrical geometry with 
azimuthal symmetry. The maximum radial extent is 
-Rmax = 20 kpc with the maximum height above the 
galactic plane z max = 5 kpc. The radiative transport 
is performed using the so-called partial intensity method 
dBaes fc Deionghdl2CM iKvlafis fc Bahcal]|ll987l) . 

Figure Q] shows the spectral energy distributions 
(SEDs) in the Galactic plane for selected galactocentric 
radii for the maximal metallicity gradient (left) and no 
metallicity gradient (right). An increased metallicity gra- 
dient reduces the UV in the inner Galaxy significantly - 
by up to a factor of 3 for A < 0.3 /im - which is redis- 
tributed into the infrared. The infrared emission for the 
inner Galaxy for the maximum metallicity gradient is a 
factor of ~ 2 higher than for the case of no metallicity 
gradient. For the outer Galaxy the ISRFs calculated for 
the two cases differ less dramatically. For the case of the 
maximal metallicity gradient the UV emission is higher 
than the no gradient case because there is less dust in 
the outer Galaxy. In turn, this results in less emission in 
the infrared than for the maximal gradient case. 

4. RESULTS 

We compare our results for the diffuse emission in the 
inner Galaxy with the new SPI data, and COMPTEL 
and EGRET data. The spectrum for each instrument is 
obtained by integrating over deconvolved skymaps. Since 
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Fig. 3. — The spectrum of the diffuse emission for 330° < I < 30°, \b\ < 15° as calculated in the optimised GALPROP model for the ISRF 
with maximal metallicity gradient, Left: with primary electrons only, Right: with secondary electrons and positrons only. Line-styles: red 
solid — 7r°-decay, green broken - IC (optical [long dash], IR [short dash], CMB [dot]), green solid — total IC, cyan solid - bremsstrahlung, 
black solid — extragalacti c 7-ray background (Strong ct al. 2004a), blue solid - total. Data points: red - EGRET and green - COMPTEL, 
as in lStrong et al J 120051) ; magenta - INTEGRAL/SPI (broke n lines: components i n fit to positronium + positron annihilation line + 
unresolved point sources; shaded region: power- law continuum) Bouchct ct al. (2008). For the SPI power-law continuum the uncertainty 
is estimated as described in the text. In this and subsequent figures, the identifier (e.g., 53J3102029RH) corresponds to the GALPROP 
version and run used; all parameters of the model are contained in the "GALDEF" parameter file for future reference and are available 
from the GALPROP website at http://galprop.stanford.edu. 



the deconvolution of the data is done based on the indi- 
vidual instrument response, cross-calibration of the data 
between the individual instruments is not an issue. 

Our comparison with the SPI data is with the dif- 
fuse e mission power- law component from [Bouchet et al.1 
(2008). We describe the procedure used to obtain the 
diffuse emission and how we estim ate the uncertainty o n 
this component. In the method of lBouchet et al.l {2008), 
a source catalogue is constructed using an iterative algo- 
rithm taking into account variable source flux contribu- 
tions using templates for the spatial morphologies of the 
interstellar emission: 8° degree Gaussian for the positron 
annihilation emission, DIRBE 4.9 /xm and CO maps for 
the continuum below and above 120 keV, respectively. 
The normalisation factor for each of these maps is ad- 
justed during the fitting procedure. Following this initial 
step, the source fluxes and template information are dis- 
carded with only the source localisations retained. The 
source position information is used in the next step of 
the analysis where the region |/| < 100° and \b\ < 30° is 
divided into cells with sizes that are chosen to optimise 
the signal-to-noise ratio per cell, while still being suffi- 
ciently small to follow the observed spatial variations. A 
likelihood fit is done using the a-priori source position in- 
formation to obtain the source fluxes and diffuse emission 
for each 'pixel' cell over the energy ranges 25-50, 50-100, 
100-200, 200-600, 600-1800, and 1800-7800 keV, respec- 
tively. This model independent "image-based" method 
establishes the extent of the diffuse emission. To ex- 
tract the diffuse spectrum with better signal-to-noise, 
the background templates (DIRBE 4.9 /im, CO) are also 



fit for each energy range. The power-law continuum is 
based on this model-dependent method but there is some 
error associated with the derived emission in this cas e 
which is not directly estimated in lBouchet et al.l {2008). 
To estimate the effect on the diffuse emission, we com- 
pare the integrated latitude profiles obtained using the 
image-based method and the fit results for the back- 
xound template maps given in Fig. 5 of lBouchet et al.1 
2008). We conclude that the intensities could be up to 
40% higher th a n the background template ones used by 
iBouchet et al.l (|2008l) to construct their spectrum. 

The primary electron and secondary positron and elec- 
tron spectra from our propagation calculations are shown 
in Figure [2j For energies < 1 GeV, the combined sec- 
ondary positron and electron flux in the ISM actually 
exceeds the primary electron flux. This is due to the 
large ratio of CR nuclei to primary electrons. The ad- 
dition of the secondary electrons and positrons increases 
by ~ 2.5 the total number capable of producing 7-rays 
via IC scattering relative to the pure primary electrons. 
Figure [3] shows the individual contributions by pri- 
mary electrons (left) and secondary electrons/positrons 
(right) to the diffuse Galactic emission. For primary elec- 
trons, the agreement with the SPI data is excellent while 
there is still some deficit when compared with COMP- 
TEL. For secondary electrons and positrons, the spec- 
trum of 7-rays is steeper below ~ 10 MeV compared to 
the primary electrons, which is a reflection of the dif- 
ferent source spectra: the primary electron source spec- 
trum is found from adjusting to the 7-ray spectrum at 
higher energies, while the secondary electron/positron 
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Fig. 4. — The total spectrum of the diffuse emission as calculated 
in the optimised model for 330° < I < 30°, |6| < 15° with the 
maximum metallicity gradient. Line and data styles as for Fig. [3] 
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Fig. 5. — The spectrum of the diffuse emission as calculated 
in the optimised model with contribution of secondary electron 
and positrons and ISRF without the metallicity gradient. Region 
330° < I < 30°, |6| < 15°. Line and data styles as in Fig. [3] 



source spectrum follows from the CR nuclei spectrum 
in the ISM. The leptonic component at low energies (IC, 
bremsstrahlung) is thus intrinsically connected with the 
higher energy hadronic component (7r°-decay). 

The components of the IC emission (Figure [3]) show 
the contributions by the ISRF components in different 
energy ranges. For primary electrons, the scattering of 
optical photons is the major contribution in the energy 
range ~ 50 MeV-100 GeV, with the infrared the ma- 
jor component below ~ 50 MeV and above 100 GeV, 
and the Cosmic Microwave Background (CMB) compa- 
rable to the infrared below ~ 500 keV. For the secondary 
electrons/positrons the scattering of the optical compo- 
nent dominates above ~ 500 keV, while the infrared is 
the major component for energies below this. Thus, 
the primary and secondary populations IC scatter the 
components of the ISRF to different hard X-ray/7-ray 
ranges. Interestingly, for secondary electrons/positrons 
the bremsstrahlung contribution is a factor ~ 2 higher for 
100 MeV to 1 GeV than the primary electron case. This 
reflects the enhancement of the electrons and positrons in 
the ISM below 1 GeV that was discussed in conjunction 
with Fig. [U 

In Figure Q] we show the diffuse emission calculated us- 
ing the optimised model with an ISRF calculated with 
a maximal metallicity gradient. Inverse Compton scat- 
tering is a major component at all energies, with ir°- 
decay more important between 100 MeV and 10 GeV, 
while the bremsstrahlung contribution is minor. The 
inclusion of secondary electrons and positrons increases 
the IC emission below ~ 100 MeV by u p to a factor 
~ 2, a n effect that was pointed out by IStrong et al.l 
(|2004rJ l. Interestingly, the agreement with the COMP- 
TEL data is improved by the inclusion of the secondary 
electrons/positrons, although the model still shows a 
deficit. Instead the SPI data are over-predicted, but 



the spectral slope is still consistent with the data be- 
low 1 MeV given the estimated uncertainties. Since the 
secondary electrons/positrons are a by-product of the 
same processes that produce the 7r°-decay 7-ray emis- 
sion at higher energies, this may indicate that the ratio 
of CR nuclei to primary electrons we use in the ISM is 
too high. A possible remedy to recover the model fit if 
the CR nuclei to primary electron ratio is reduced could 
be that we have simply underestimated the optical com- 
ponent of the ISRF. If the CR nuclei flux is reduced to 
improve the agreement with SPI the model emission in 
the EGRET energy range would also be reduced. The 
emission in the MeV energy range comes from secondary 
electrons/positrons IC scattering the optical component 
of the ISRF, and primary electrons IC scatter the same 
component to GeV energies. Increasing the optical ISRF 
would simultaneously increase the emission in the MeV 
and GeV range, recovering the agreement of the model 
with EGRET data and possibly further improving the 
agreement with COMPTEL. 

To test the dependence of the sub-MeV emission on the 
assumed ISRF model, we calculate the diffuse emission 
for the case of no metallicity gradient. Figure[5]shows the 
spectrum of the diffuse emission for this case. The major 
change in the ISRF is a ^30% reduction in the infrared 
emission which is due to the smaller amount of dust in 
the inner Galaxy when using this ISRF model. This re- 
sults in a drop of ~10— 15% for the predicted intensity 
of the diffuse emission below ~1 MeV than the maxi- 
mal gradient case. The optical emission increases only 
slightly, while the CMB emission stays the same. The 
total spectrum of the diffuse emission does not change 
significantly under this variation of the ISRF showing 
the robustness of the calculations. 

We point out that the reduction in the infrared signif- 
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Fig. 6. — The spectrum of the diffuse emission as calculated in the optimised model with contribution of secondary electron and positrons 
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icantly reduces the contribution by secondary electrons 
and positrons to the total emission below a few hundred 
keV. In this energy range the emission by primary elec- 
trons is also reduced, but by a smaller amount because 
there is still a contribution by IC scattering of the CMB 
(see Fig. [3|). Since the CMB is known, and if the in- 
frared component of the ISRF is low, the emission below 
a hundred keV traces the primary electron spectrum in 
the ISM. 

We have made our principal comparison with the spec- 
trum for \b\ < 15° since this is the nominal ra nge for the 
spectrum presented in iBouchet et al.l (|2008). To illus- 
trate the effect of different latitude ranges, we calculate 
the diffuse emission for \b\ < 10° and \b\ < 5° and com- 
pare with the data for these ranges. The results for these 
latitude ranges are shown in Fig. [5] Since the extraction 
of the diffuse component from the SPI data relies upon 
templates whose latitude distribution may not exactly 
match the true distribution of the emission some of the 
signal may be absorb ed into the baseline (see Fig. 5 of 
IBouchet et al.l (|2008f )). This introduces further uncer- 
tainty when comparing reduced latitude ranges which is 
difficult to quantify better than we have already done. 
The model emission is qualitatively similar to the data 
for the reduced latitude ranges, which points the way to 
using the IC emission as a template in future analyses of 
the SPI data. 

5. DISCUSSION AND CONCLUSIONS 

The agreement of the single model over the whole en- 
ergy range from SPI data at low energies, to the EGRET 
data at high energies, is rem arkable. Note that we are 
using the "optimised model" ([Strong et al . 2004b) which 
has higher primary CR electron fluxes than observed lo- 
cally, and also higher CR nuclei fluxes, as required to 
fit the 7-ray data above 30 MeV. With this model we 



predict too much emission below 1 MeV but seem to 
reproduce the spectral slope. Reducing the CR nuclei 
source spectrum to improve the sub-MeV agreement is 
a possibility. However it cannot be reduced too much 
since the local CR antiproton fluxes must still be repro- 
duced by the model 4 . Increasing the optical component 
of the ISRF could improve the agreement at MeV and 
GeV energies. Variation of the primary electron injection 
index below a few GeV is another possible remedy to the 
over-production of diffuse emission below 1 MeV. This 
requires a different source spectrum than is presently 
used in the optimised model to sufficiently reduce the 
sub-MeV diffuse emission in order to be consistent with 
the SP I spectrum. If inste ad we use the "conventional" 
model ([Strong et al.ll2004bT ) the situation will not be im- 
proved: the sub-MeV emission will be lower, but the 
agreement with the COMPTEL and EGRET data will 
be substantially worse. 

There is still room for a contribution from popula- 
tions of unresolved compact sources, particularly anoma- 
lous X-ray pulsars and radio pulsars, which may have 
very hard spectra extending to a few hundred keV 
(jKuiper et aLl feOOB). They may be responsible for the 
apparent peak near b — 0° in the IBouchet et al.l (|2008[ ) 
latitude profiles. 

The hard X-ray continuum is consistent with the pre- 
dictions in both intensity and spectral index. How- 
ever, the uncertainties in the model are still consider- 
able: the distribution of CR sources and gas in the in- 
ner Galaxy which affect both the primary and secondary 



4 The BESS- Polar flight of 2004 IIHams et alJ2007l) revealed that 
the CR antiproton flux is somewhat lower than previous measure- 
ments with lower statistics. If these data are confirmed, the CR 
proton spectrum will become better constrained. In turn, it will 
lead to a somewhat smaller flux of secondary leptons and improved 
agreement with the INTEGRAL data. 
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electrons/positrons, and the optical and infrared part of 
the ISRF (the CMB is of course known exactly). In fact, 
our optimised model overpredicts the SPI data, which 
could simply reflect these uncertainties, but the agree- 
ment in the spectral shape gives confidence that the 
mechanism is correctly identified. 

There is still an excess in the COMPTEL energy range 
between 1-30 MeV which is to be explained. The con- 
tribution of the positron anni hilation in flight may con - 
tribute in this energy range (jBeacom fc Yuksel |2006) , 
but it has to be tested against the intensity of the 511 
keV line and positronium continuum. 

From our modelling, we find that the total rate of sec- 
ondary positron production by CRs in the whole Galaxy 
is ~ 2 x 10 42 s _1 in the optimised model. The conven- 
tional model gives a factor of —2 less positrons. These 
values are ~ 10% of the positron annihilation rate ~ 
1.8 x 10 43 s _1 as derived from I NTEGRAL observation s 
of the 511 keV line emission (jKnodlseder et all 12005) . 
The current CR flux of positrons is not sufficient to ac- 
count for the observed annihilation rate. A CR origin for 
the 511 keV annihilation line can still be reconciled with 
the production rate if CR intensities in the past were 
higher. 

Our work illustrates the intrinsic connection between 
the diffuse Galactic 7-ray emission in different energy 
ranges. Inverse Compton emission by CR electrons and 
positrons on starlight and infrared radiation are the most 
important components of the hard X-ray and 7-ray emis- 



sion in the 100 keV to few MeV range. A considerable 
proportion of this emission is produced by secondary 
electrons and positrons, the spectrum of which depends 
on the CR nuclei spectrum at energies ~ few GeV and 
higher. These CRs also produce 7r°-decay 7-rays that 
dominate the emission in the GLAST range from 100 
MeV to -10 GeV. Hence, GLAST observations of the 
7r°-decay diffuse emission will also constrain in the fu- 
ture the contribution by secondary electrons/positrons to 
the diffuse 7-ray emission in the SPI energy range. With 
the secondary electrons/positrons fixed, SPI observations 
probe the IC emission of primary CR electrons with en- 
ergies < 10 GeV scattering the infrared component of 
the ISRF and the CMB. This will provide information 
on the low energy spectrum of primary CR electrons and 
the infrared component of the ISRF. In turn, since most 
of the diffuse 7-ray emission between ~10 GeV - 10 TeV 
is produced via IC scattering of primary electrons on 
the same starlight and infrared photons, this provides a 
connection to observations of diffuse emission at TeV en- 
ergies by HESS (lAharonian et al.l l2006t) and MILAGRO 
(lAbdo et al.ll2007lK 
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